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ABSTRACT 


Second-order steady wave forces on a restrained, submerged, 
circular cylinder are studied experimentally as a function of 
cylinder depth and wave frequency, The experimental results are 
compared to forces calculated from a linearized potential theory, 
It is found that in all cases the experimentally observed forces 
are higher - in some cases by an order of magnitude - and are not 


attenuated with depth as rapidly as predicted by the theory. 
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1, INTRODUCTION 


Since very early in written history warriors have 
attempted to escape the sea’s surface and OREEE submarine 
warfare. Herodotus (460 B.C.), Aristotle (332 B.C.), and 
Pliny, the elder (77 A.D.) mention determined attempts to 
build submersibles, Alexander ae Great (356 to 323 B.C.) 
is the first person know to have descended into the sea in 


a vessel of any eds 


Since that time submarines have been invented which 
were propelled by oars, sails, treadles, hand-operated screws, 
clockwork, springs, steam stored in tubes, chemical engines, 
compressed air, stored gases, electric motors, and nuclear 
power, Over three hundred years azo Mother Shipton, a famous 
English prophetess, predicted the coming of the submarine when 
writing, "underwater man shall walk, shall ride, shall sleep, 
saaii talk,” Her mrennhesy has come true, but in carrying 
out her prophesy man has discovered that even though he may 
be an underwater man in many senses, he has not completely 


escaped the forces of surface waves. 


Man was early to recognize the forces of the sea to be 
powerful, at times destructive, and inevitable. He set out to 
study them, attempt to predict them, and to design vessels 
which would withstand and take advantage of them. Studies of 


forces on vessels penetrating a free surface are too numerous 
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to mention here, 


The United States Navy commissioned the USS Holland, its 
first submarine on April 11, 1900. The military importance of 
submarines has increased at a staggering rate since that time. 
Nore recently, commercial and scientific interests have awakened 
to the vast resources uncer the sea, hence the use of non- 
military submersibles has begun to increase at a rate that may 
exceed that of its military cousin. As more and more use is 
made of undersea vessels, they become larger -- the largest 
knotm maneuverable undersea vessel to date is about ten meters 
in diameter, greater than one hundred meters in length, and 
displaces in excess of six thousand tons -- and yet they must 


be precisely controlled, 


In order to be able to provide this precise control, the 
forces acting on these vessels must be knotn, One source of 
such forces is surface waves. A number of studies have been 
made of the ate induced forces on underwater bodies, but most 
of these studies have been concerned with an underwater body 
moving with some translational velocity with respect to the 
fluid. Examples of these studies are some of Havelock's 


3 


papers, Kim's ~ paper, and eeierecenaee paper. 


The present work will concern itself with a body that is 
not moving, but which is rigidly restrained. Kulin? studied the 
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the effect of a single intumescence which subjects. the fluid to 
acceleration followed by deceleration. He studied (1) the 

extent local acceleration causes the drag coefficient to depart 
from steady state values, (2) the extent that friction affects 
potential flow values of inertial coefficient, and (3) the extent 
of force dependence upon preceding flow history, as well as on 
instantaneous conditions. He found that the forces could not be 
represented by a single universal total resistance coefficient, 
and that the ratio of the distance of fluid motion to object 

Size Was a significant parameter related to vortex formation 


during fluid deceleration, 


lita aisle: Ogilvie? applied and extended the form of 
solution that had been developed me Wetec! in 1950 to calculate 
the first-order oscillatory force and the second~order steady 
force in the following situations: (a) a submerged cylinder is 
restrained from moving under the effect of incident sinusoidal 
waves; (b) a submerged cylinder is forced to oscillate sinusoidally, 
in otherwise calm water; (c) a cylinder, which is neutrally 
buoyant, is allowed to respond to first-order oscillatory forces, 
_He proves that subject to his assumptions a knowledge of the first- 
order potential supplies information sufficient to solve these 
problems, This paper deals with an attempt to experimentally 


verify his results for case (a), the restained cylinder. 


Assume that a circular cylinder is located under a free surface, 





(see figure 1.) such that the surface of the cylinder is 


specified by 
S(x,y) = x* + dy + ne ae toe (1) 


The undisturbed free surface is taken as the x-axis and 


the instantaneous free surface will be specified by 
y - Y(x,t) = 0. (2) 
The required velocity potential, O(x,y,t), satisfies 


D, + Dy, =0 for y < Y(x,t), (3) 


ance Gs) > 0, 


io = 1 +, = 0 Ore = Y(X,t), (4) 
gy + De + 4( mf + D ) = (0) on y = NE Crop oy (5) 
and 


(vo ye sytoo. = © on S(x,y) = 0, (6) 


- 


t 
plus appropriate conditions at infinity. The time average, Z(t), 


of the second-order steady force can be shown to be 
¢ ae t 
gO) = J Ge) = aA) 
T 
es olay Ty | de (7) 
= -zlap je LJP ae y ee 7 
“7 


or if we let 0(x,7,t) be the real part of a function of 





a complex variable, f(z,t), where z = x + iy, then 


; Tr 
RL) = “Hap [eo Ff (2,t)f'( 0 |,_,29 (8) 
r 
= ~i2mp2(HO) e ~*" Q(ka,kh). (8°) 


k is the wave number defined as 


k =o4/¢, (9) 


and HO is the amplitude of the incident waves, 


Define steady vertical force coefficient,f, as 


b-PR. | (10) 


2mpe(SOy 


From (8') and (10), it is seen that 
B =e" Q(ka,kh), (11) 


It is this quantity obtained by Ogilvie that will be 


compared with experimental results in Section IV. 





Figure 1, Geometry of the problem 
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II, TEST APPARATUS AND PROCEDURE 
IT.A. TEST APPARATUS DESCRIPTION 


This experiment was conducted in the wave tank in the M.I.T. 
Department of Naval Architecture and Marine Engineering Hydro- 
dynamics Laboratory. The working section of the tank is ten 
feet long, one foot wide, and six inches deep. The wave generator 
is a vertically oriented pusher, hinged at the bottom, and driven 
by a push--rod connected off center on a flywheel, which is in turn 
driven by a variable speed electric motor. The observation section 
was located about equidistant from the two ends, An energy 
absorbing beach made of rubberized horse hair was located at 
the end opposite the wave generator, This beach was not completely 
successful in preventing reflected waves, This will be considered 


maeoeetion 11.8, 


Wave height was measured by a capacitance probe mounted on 

a carriage wnich was driven by a variable speed electric motor 
and ran on tracks along the top of the tank, The output of the 
probe was fed through a buffer amplifier to a Hewlett Packard 
Model 7702 B recorder, The carriage was moved at a constant 
speed before each force observation to obtain data for cal- 
culating the reflected to direct wave ratio and then placed 

at a predesignated point as close as possible to the cylinder 


to measure wave height during the period of force measurement. 
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The cylinder itself, a one inch diameter lucite rod 11 15/16 
inches long, was suspended below the surface by a 3/32 inch 
diameter rod tapped into the cylinder at points 1/8 inch from 
each cylinder end, Any effect on vertical forces by these thin 
supporting rods was assumed to be negligible. The 1/32 inch 
clearance between each end of the rod and the tank wall is 
considered to be smaller nora any wall boundary layer thickness, 
and it is felt that there was no appreciable departure from the 
two dimensional condition. The tops of the supporting rods 
were Pesadea and attached to a yoke which was restrained so 


that its only degree of freedom was in the vertical direction, 


Vertical force on the yoke was measured by two Schaevitz 
Engineering Wedel No. FTA 3 Force Transducers, one on each end, 
The electrical outputs of these transducers were connected and 
processed as shown schematically in figure 2. As seen from the 
figure the outnut voltage was directly proportional to the time 


integral of the steady force component, a. 
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Ii.B. PROCHDURE 


As mentioned in Section II.A., the energy absorbing beach 
was not completely successful in preventing reflected waves and 
it was considered necessary to take these reflected waves into 
account. A method of doing this by the use of a standing wave 
ratio, p, is developed in Appendix A. That result is given in 


(A7) and is.repeated here. 
p ~ (oe (12) 
(1 + p? 


The values of HX and HM used in (A8) were obtained by 
allowing the capacitance probe carriage to move at constant 
speed along the tracks as mentioned in Section II.A. and observ- 


ing the maximum and minimum values of wave height recorded, 


The wave height recorder was calibrated by changing the 
still water level by a known amount and observing the change 
in recorder level. Because of the sensitivity of a capacitance 
probe, the calibration was repeated frequently and the appropriate 


calibration factor was used in each calculation of B. 


Force measurement calibration was accomplished by first 
opening the shorting switch showm bypassing the integrating 
amplifier in figure 2, for a known period of time, TD, and 
observing the intergrated drift voltage, ED. The switch was then 


evocee send 


the evlinder sunnorting 





yoke, The switch was reopened for a known period of time, TS, 
and the integrated signal voltage, ES, was observed, The cali- 


bration factor, AE, was then determined by the following formula: 
AE = Y/E ee (13) 

where 
eeen(no/tS) = (2D/TD). (14) 

This factor was found to remain constant at 0,036 pounds per volt. 


: t 
The unknowm steady wave force, Y(t), was found by the 
following equation, using values of ES, TS, ED, TD, HM, and HX 
observed for various cylinder depths, h, and incident wave 


frequencies, F, 


ei ig 
YALC) = (E) (AB) (15) 
BOwas calculated using (10) and (15) and was calculated by (12). 


Calculations were performed at the M.I.T. Computer Center 


using the program in Appendix B. 
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ITI. RESULTS 


Experimental data is tabulated in Appendix C. Values 
of #calculated from this data using (12) are tabulated on 


the following page. 


These values of @ are presented graphically in fisures 
3. through 9. along with curves for ka = 0.2 and 0.5 from 


Ogilvie's theory. 


In drawing the experimental # curves, curves for the next 
higher and lower values of ka were considered as well as @ 
values for the curve being dratm,. As a result the experimental 


curves for some ka values do not pass through all data points, 


The individual @ curves from figures 3. through 9. are 


reproduced in figure 10. for ease of comparison. 
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Figure 3. Steady vertical force coefficient, ka = 0725 
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Figure 4, Steady verticai force coefficient, ka = 0. 
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Figure 5. Steady vertical force coefficient, ka 
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Figure 6, Steady vertical force coefficient, ka = 0.4 
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Figure 7, Steady vertical force coefficient, ka 
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Figure 8, Steady vertical force coefficient, ka = 0.5 
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Figure 9, Steady vertical force coefficient, ka = 0.55 
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Figure 10. Steady vertical force coefficients 
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_LiV. DISCUSSION OF RESULTS 


Reference to figures 3. through 9. shows that the 
experimental curves of steady vertical force coefficient are 
in all cases higher -- in some cases more than an order of 
magnitude ~~ than the theoretical curves, and in all cases 
eeepu tor ka = 0,25 and 0.3 the slope is less steep than the 


theoretical curves, 


It is also noted that the experimental curves all have 
inflection points between the third and fourth data points -- 
at an h/a ratio of about 4, From figure 10. it is easily 
seen that this inflection is much more pronounced in the curves 
for smaller values of ka -- longer wave lengths. Since for the 
longest wavelength investigated the tank depth was over a half 
Wave length it was expected that it would have an infinite 
Eeaeerive deptn; however, a finite depth effect seems to be a 


plausible explanation for these curve inflections, 


ieee 





Ve. CONCLUSIONS AND RECOMMENDATIONS 


From Sections III and IV it may be concluded that Ogilvie's 
potential theory does predict the general behavior of the steady 
vertical force coefficient, but that in all cases it predicts 
2 force smaller than that measured experimentally and that the 
forces will be attenuated more rapidly with depth than that 


experimentally observed. 


Because of the curve inflections noted in Section IV and 
Since smaller values of ka are more likely to be encountered at 
sea, 1.€., a one hundred meter wave length and a ten meter 
diameter cylinder give a ka of 0.314 -- most waves of interest 
are longer and most diameters are smaller, hence most ka‘s of 
interest are smaller, it is felt that an extension of this 
experiment to smaller values of ka in facility with an infinite 


effective depth and no reflected waves would be rewarding. 
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APPENDIX A 
Theory Modification to Account for Reflected Waves 


Using the ‘circle theorem" of oneethonesone the complex 


potential function, f(z,t), in (8) can be represented by 


Hot +kz) kh ee | 
jcZ.t) = A ee Zz + ih (A1) 


where 


A = (HO) elo. (A2) 


Using (A1), (8) and (10), is found to be 
-2kh 

B=e (ka) I, (2ka). (43) 

The complex potential function for a reflected wave 
traveling in the opposite direction with an unknotm phase 
angle, Y, and a mamitude of p times the direct wave can be 
given by 

 a¥[ -a@t - kz) -kh se - ke 

f.(z,t) = phe |e +e oe Z- ih)jl, ( ALL) 

Now the complex potential function for a direct and reflected 


wave existing simultaneously is given by 
tet f (z,t) + f,(z,t). (A5) 
Using (A5), (8) and (10), BO is found to be 


~2kh 
Bor Glare Je)” 0 (ka) I, (2ka). (A6) 





Consequently 


B = B/G + Pp) (47) 


where B would be the steady vertical force coefficient if 
there were no reflected wave present, 80 is the steady vertical 
force coefficient with a reflected wave present, p is the 
standing wave ratio given o 

p= HX - HM (A8) 

EDS ce au 

and HX and HM are the maximum and minimum wave amplitudes 
that exist as a result of the summation of the direct and 
reflected waves. HO used in (A2) is the average wave amplitude 
given by 

HO = HX + HM, (A9) 

Ze 

The complex potential functions given by (Al) and (A4) do 
not satisfy the free surface conditions (4) and (5), hence 
(A3) and (A65) can be expected to give reasonable results only 
if ka<e< 1 and/or 2kh »> 1. However, it is assumed that (A7) is 


valid throughout the interest range of ka and kh values, 


BOGS. 





APPENDIX B 


Computer Program for Steady Vertical Force Coefficient Calculation 


Piclor 
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APPENDIX C 


Experimental Data Tabulation 
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